Elastic instability in stratified core annular flow 
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We study experimentally the interfacial instability between a layer of dilute polymer solution and 
water flowing in a thin capillary. The use of microfluidic devices allows us to observe and quantify in 
great detail the features of the flow. At low velocities, the flow takes the form of a straight jet, while 
at high velocities, steady or advected wavy jets are produced. We demonstrate that the transition 
between these flow regimes is purely elastic - it is caused by viscoelasticity of the polymer solution 
only. The linear stability analysis of the flow in the short-wave approximation captures quantitatively 
the flow diagram. Suprisingly, unstable flows are observed for strong velocities, whereas convected 
flows are observed for low velocities. We demonstrate that this instability can be used to measure 
rheological properties of dilute polymer solutions that are difficult to assess otherwise. 

PACS numbers: 83.80.Rs, 83.60.Wc, 83.85.Cg 
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Polymer solutions exhibit purely elastic flow instabili- 
ties even in the absence of inertia [l| . The almost ubiq- 
uitous ingredient of such an elastic instability is the cur- 
vature of streamlines: polymers that have been extended 
along curved streamlines are taken by fluctuations across 
shear rate gradient in the unperturbed state which, in 
turn, couples the hoop stresses acting along the curved 
streamlines to the radial and axial flows and amplifies 
the perturbation @, Q . Flat interfaces between two fluids 
with different viscoelastic properties can also become un- 
stable 0Bl due to the normal stress imbalance across the 
interface. These instabilities often occur in coextrusion 
where different polymers are melted in separate screw ex- 
truders and then flown simultaneously in the extrusion 
nozzle. Undesirable wavy interfaces are sometimes ob- 
served between the adjacent polymer layers both during 
the flow and in the final product Q. Since these insta- 
bilities set severe limits to the industrial processes such 
as film or fiber fabrication, they have been extensively 
studied before 0, H| ■ Althoughprevious experiments and 
theory agree reasonably well 0, Q , a comprehensive de- 
scription of the flow is still lacking fioj . 

In this Letter we propose a quantitative explanation 
for various flow patterns observed in purely elastic inter- 
facial instabilities. We perform a set of original exper- 
iments on co-flow of a polymer solution and water and 
map the full flow diagram. In contrast to previous exper- 
imental studies dealing with macroscopic flows of molten 
polymers, we focus here on flows of dilute polymer solu- 
tions and water in microfluidic devices. The advantage of 
using dilute polymer solutions is that their elastic prop- 
erties can easily be tuned by dilution and that most of 
the theoretical work has been performed for models rep- 
resenting dilute polymeric fluids. The use of microfluidic 
flow geometries offers simple control and visualisation of 
the flow. 

We observe that above some flow rate the interface 



between the polymer solution and water becomes wavy. 
Surprisingly, for relatively low velocities, the instability is 
convected downstream, while a stationary unstable flow 
is observed at high velocities. This behaviour is due to 
the interplay between advcction by the mean flow and 
the growth of a perturbation and is in close agreement 
with the model we propose. Our work opens the route to 
use such an elastic instability in straight channels to pro- 
mote mixing in microfluidic devices . We demonstrate 
that it is also a new way to measure rheological proper- 
ties of weakly elastic polymer solutions where very few 
techniques are available. 

Our experiments are performed in a microfluidic de- 
vice made of nested glass capillaries: an inner capil- 
lary of square cross-section with tapered nozzle (section 
~ 300 urn) nearly perfectly fits into a cylindrical capil- 
lary (R c — 400 urn) that carries the outer fluid; it offers a 
simple way to self-center and align the capillaries (l2l.[l3j]. 
The length between the nozzle and the outlet of the de- 
vice is set to L = 6 cm. The two co-flowing fluids are 
injected with precision syringe pumps at flow rates Qi 
and Q e for the internal and external rates respectively. 
The microfluidic chip is positioned vertically on a stan- 
dard microscope (mounted accordingly) in order to pre- 
vent the effect of gravity. The observations are carried 
out with a fast camera (Miro Phantom). Our working 
fluid is a solution of Poly(VinylAlcohol) (PVA) of mo- 
lar mass M w = 196000 g/mol in water at concentrations 
of 3.25, 5, 6 and 7.5 % wt/wt (for which c > c* w 1% 
wt/wt). The measured values of the shear viscosity of 
the solutions are given in Table U and do not depend on 
the shear rate up to 10 2 s _1 . The viscosity of water is 
taken to be % ~ 10~ 3 Pa - s at 20 °C. 

We observe flow patterns that depend strongly on the 
flow rates. Fig. Q] shows a typical flow diagram in the 
(Qi, Qe) plane for water and the semi-dilute PVA solu- 
tion as the inner and outer fluids, respectively. For low 
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flow rates, we find straight jets (+) that extend up to 
the outlet of the device, while at higher flow rates the 
interface between the two fluids becomes unstable. At 
intermediate flow rates, jets are straight from the inlet 
up to some distance downstream where varicose undula- 
tions set in. We call them advected wavy jets (o) and note 
that this distance decreases with the flow rates. At yet 
higher flow rates, jets are wavy through the whole set- 
up (•). We note here that the advected wavy and wavy 
jets are, most probably, the same dynamical state, the 
only difference is that the latter sets in before and for- 
mer after some arbitrary lengthscale L c . We distinguish 
between the two states in order to be able to measure 
the relaxation time of the polymer solution, as will be- 
come apparent later. Unless mentioned otherwise, we set 
L c = 10R C . 

In order to identify the origin of the instability, we 
first calculate the laminar flow profile in the system. We 
neglect both inertial effects and molecular diffusion pro- 
cesses since the Reynolds number is small (Re ~ 0.1) 
and the Peclet number is large (Pe > 10 4 ) in our exper- 
iments. The two miscible fluids thus flow side by side 
without mixing and exhibit a constant effective surface 
The motion of water is then described by 
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tension 

the Stokes equation while the polymer solution obeys the 
Oldroyd-B model [H: 



-Vp + rj s Av + V-S = 0, 



(1) 



where £= <9 t £ + w- V£-(W) T ■ S-E ■ Wis the upper- 
convected derivative [l5j |. Here, v is the velocity and p is 
the pressure in the fluid, S is the polymer contribution 
to the stress tensor; ry s is the viscosity of the solvent, and 
t and rj p are the Maxwell relaxation time of the polymer 
and the increase of viscosity due to the polymer chains 
respectively. The total shear viscosity of the polymeric 
solution is thus r\ e = r) p +i~i a . We enforce no-slip boundary 
conditions at the solid-liquid interface and the continuity 
of the velocity and of the tangential stress at the inter- 
face between the two fluids. In the unidirectional laminar 
flow of Fig. [TJ the pressure drop AP between the nozzle 
and the outlet of the capillary is the same in both flu- 
ids and is related by Eqs. ^ to the flow rates and the 
relative position of the interface between the two fluids 
x = Ri/R c : 



AP 



nRi(l 



(2) 



1 + nijy-, and m = r] e /r)i is 



where x — . 

V a—l+m ' V Qc 

the viscosity ratio. 

In the left part of Fig. [2] we redraw our experimen- 
tal data in the (x, AP/ L) plane and observe that stable 
flows occur at low while unstable flows occur at high 
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FIG. 1. (Left) Experimental flow diagram of the 5% PVA 
polymer solution (external fluid) and dyed water (internal 
fluid) as a function of the respective flow rates. (Right) Flow 
patterns observed in the microfluidic chip: stable straight jets 
(+), wavy jets (•) and advected wavy jets (o). 




FIG. 2. (Left) Flow diagram in the (x,AP/L) plane [same 
data as in Fig. [T] stable straight jets (+), wavy jets (•), ad- 
vected wavy jets(o)]. (Right) Same data compared with the 
kinetic criterion Q: jets with a straight part longer (x) and 
shorter (•) than L c = 4.8 mm. The solid line is calculated 
from with the polymer relaxation time r = 9 ms. The 
two dotted lines with r = 7 ms and r = 11ms bracket the 
uncertainty. 



pressure drops. This observation allows us to dismiss the 
Rayleigh-Plateau mechanism as a possible origin of the 
instability that would be triggered by the effective surface 
tension between the two miscible fluids 
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this would be contradictory to recent experimental and 
theoretical results [l3| that explicitly demonstrate that 
droplets and wavy jets (absolutely unstable states) occur 
at low pressure drops while straight jets dominate at high 
pressure drops. 

We also exclude the viscosity stratification of the flow 
as an origin of the instability EH • We have repeated 
the same experiment with a glycerin solution with r] e = 
0.1 Pa - s and we have never observed unstable inf erf aces. 
Since the viscosity-contrast based instability is of inertial 
origin, it should develop at higher Reynolds numbers. 
However, it does not play a role at our flow conditions. 

We are thus left with a purely elastic instability 
driven by the constrast of the normal stresses across 
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FIG. 3. State diagrams for the 3.25 %wt PVA polymer solu- 
tion determined for several reference distances L c . A single 
value t = 8 ms is used to fit all the state boundaries. 

the interface This instability was studied analyt- 
ically by Chen Q and by Chen and Joseph [l9[ as- 
suming low Reynolds and high Peclet numbers as in 
Eq. ©. They performed the linear stability analysis of 
the flow with respect to small axisymmetric perturba- 
tions oc exp (ikz + tot) with lo being the complex growth 
rate and k setting the wavelength of the perturbation. 
They found that long-wavelength perturbations are al- 
ways stable In the opposite limit k — > oo, the flow is 
always unstable [l9| and for weak elasticity of the solu- 
tion, the dispersion relation can be approximated by: 

m(m — 1) / AP xR c \ 2 , AP R 2 r , 9 , 

oj = — t - ik -(1 - x ). 

(1 + m) 2 \ L 2 Ve J L 4 Ve K ' 

(3) 

An important feature of this dispersion relation is that 
the real part of the growth rate is independent of the 
wavelength. Full numerical linear stability analysis (to be 
published elsewhere) confirms that the dispersion curve 
is practically flat and positive only becoming negative for 
very small k's. This implies that almost all wavelengths 
become unstable with an identical growth rate and the 
question of the wavelength selection cannot be answered 
based on the linear theory. Intriguingly, the use of a 
more complex rheological model would permit to raise 
this degeneracy 0- 

In order to describe the convective nature of the insta- 
bility we propose a simple kinetic criterion that captures 
most of our experimental observations. We assume that 
the instability sets in very close to the nozzle and is grow- 
ing on the typical timescale r r = 1/Re(uj) while being 
advectcd downstream with the velocity of the interface 
U. The typical development length of the instability is 
then L = Ur r , and the boundary between the advected 
wavy and wavy jets is given by L = L c . In terms of the 
applied experimental parameters, this criterion reads: 

_ AP L c t (m - l)m _ 1 - x 2 
* = ~T~ (1 + m) 2 = x 2 ' () 

This criterion offers a few interesting predictions. First, 
the dependence upon the pressure drop is quite surpris- 




FIG. 4. Comparison between the experimental state diagram 
and the stability criterion ([4]) (solid line) for several polymer 
solutions (different concentrations). Solid symbols (•) corre- 
spond to unstable flows and crosses ( x ) to stable jets on the 
length L c . 

ing. Eq. Q implies that for a given ratio of the flow 
rates (which is independent of AP), the advected wavy 
jets (L > L c ) occur at low pressure drops, while the wavy 
jets (L < L c ) should occur at high pressure drops, as ob- 
served experimentally. Since the velocity of the interface 
scales linearly with the pressure drop while the destabilis- 
ing forces due to the viscoelastic normal stresses scale as 
AP 2 [lij], the instability moves closer to the inlet upon 
increase of the pressure drop. This is in contrast with 
the Rayleigh-Plateau instability due to the surface ten- 
sion for which the opposite order of dynamical states is 
observed Moreover, the higher the pressure drop, 

the smaller the polymer relaxation time must be to reach 
the condition L = L c , that turns out to be a fruitful way 
to measure t, as we show below. We also note that this 
is a purely elastic instability as it vanishes in the Newto- 
nian limit r = 0. Finally, we observe that the criterion 
(|4]) is indepedent of the size of the capillary that sug- 
gests that this instability will also exist in nanofluidic or 
macrofluidic devices provided the incrtial forces are kept 
small. 

We now compare this kinetic criterion to our experi- 
ments. The only unknown quantity in @ is the polymer 
relaxation time r and we first use it as a fitting parame- 
ter. The right panel of Fig.[2]shows that the theory based 
on the single-relaxation time Oldroyd-B model agrees 
reasonably well with the experiments. For the 5% so- 
lution we have extracted r = 9 ± 2 ms. Some discrepancy 
however is apparent at small radii of the Newtonian core 
x. One possible source of the discrepancy is the approx- 
imate nature of the dispersion relation (j3|). It is derived 
for short- wavelength perturbations while we use it for dis- 
turbances of intermediate wavelength. This only makes 
sense if the actual dispersion relation is flat for most of 
the fc's. The full numerical linear stability analysis shows 
that for small x the dispersion relation is less flat than 
for large values of x which possibly explains the discrep- 
ancy between theory and experiments in Fig. [2j Another 
possibility is that at small x we observe a decrease of the 
relaxation time with the local shear rate (small x corre- 



4 



TABLE I. Comparison of the polymer relaxation time ob- 
tained from the instability study (r), the drop detachment 
experiments (r a ) and shear rheometry done at 7 = 10s _1 (t w ) 
for several concentrations of PVA (M w = 196000 g/mol). 



wt% 


r] e (Pa-s) 


t (ms) 


T a (ms) 


t w (ms,7 = 10s 1 ) 


7.5 


0.65 


120±20 


70± 20 


100 


6 


0.25 


60 ±10 


20± 10 




5 


0.1 


9 ±2 


15± 10 


12 


3.25 


0.04 


8 ±2 







spond to large pressure drops and thus high shear rates 
at the interface). 

Next we show that the model is actually self-consistent 
as the measurement of r is quite insensitive to the choice 
of L c . In Fig. [3] we plot the stability diagram of the 
same system with various reference distances L c . A sin- 
gle value of the polymer relaxation time is able to reason- 
ably fit most of the data. Once again, the discrepancies 
at small x are probably due to the reasons mentioned 
above. 

The experimental results for several polymer solutions 
are summarised on the master phase diagram in Fig. [4] 
Clearly, the simple criterion ((4]) is remarkably successful 
in predicting the transition between advected wavy and 
wavy jets. In Table|T]we provide the values of t extracted 
for several concentrations of the polymer. To check these 
values, independent measurements of the polymer relax- 
ation time have been performed. First, we use a version 
of extensional rheometry described by Amarouchene et 
al. . When a drop of a polymer solution detaches from 
a capillary tube, a long-lived cylindrical neck is formed. 
Initially, it thins according to a power law for all liquids, 
and then further, exponentially in time if the liquid is 
viscoelastic. The decay time of the exponential thinning 
is directly proportional to the characteristic time of the 
polymer r a . We have also performed conventional rhe- 
ological measurements of the polymer relaxation time. 
The data are quite noisy since the solutions are not very 
viscous and only weakly elastic. In Table U we report 
the values t w measured at the shear rate of 10 s _1 which 
is similar to the shear rates in our coflow experiments. 
These values are in a reasonable agreement with the val- 
ues r extracted from the onset of wavy stationary jets. 
Conventional rheometry also shows that the polymer re- 
laxation time t w is a decreasing function of the shear rate 
providing support for our explanation of the discrepan- 
cies in Figs. [5] and [3] 

In this Letter, we have demonstrated that the purely 
elastic interfacial instability that exists between two flu- 
ids with different viscoelastic properties can not only 



be explained on a kinetic basis but also used to mea- 
sure the polymer relaxation time in weakly elastic poly- 
mer solutions which is quite difficult to obtain otherwise. 
The technique also has the potential to measure shear- 
thinning of the solution at high shear rates. We have 
performed the fitting procedure described above indepen- 
dently for every x and thus obtained the curve r vs shear 
rate. The results look promising and presently we are 
working on the extension of the technique to more com- 
plicated constitutive relations than the Oldroyd-B model 
that will allow us to extract full non-linear rheology even 
for very dilute polymer solutions. 

We would like to thank Dr. Christian Wagner and 
Christof Schafer, Saarland University, for their help with 
the rheometry. AM acknowledges support from the 
Royal Society of Edinburgh/BP Trust Personal Research 
Fellowship and EPSRC Career Acceleration Fellowship 
(grant reference number EP/I004262/1). LOF acknowl- 
edges support from the Aquitaine Council and from ANR 
Pnano project MicRheo. 



[1] 

[2] 
[3] 

K 
[5 

[6 

It 

[8 

[9 
[10 

[11 
[12 

[13 

[14 
[15 

[16 
[17 
[18 

[19 

[20 



Agassant, Int. 
Agassant, Int. 



annie.colin-exterieur@eu.rhodia.com 

R. G. Larson, Rheo. Acta 31, 213 (1992). 

E. S. G. Shaqfeh, Annu. Rev. Fluid Mech. 28, 129 (1996). 

G. McKinley, P. Pakdel, and A. Oztekin, J. Non-Newton. 
Fluid Mech. 67, 19 (1996). 

K. P. Chen, J. Non-Newton. Fluid Mech. 40, 155 (1991). 
E. J. Hinch, O. J. Harris, and J. M. Rallison, J. Non- 
Newton. Fluid Mech. 43, 311 (1992). 

H. J. Wilson and J. M. Rallison, J. Non-Newton. Fluid 
Mech. 72, 237 (1997). 
R. Valette, R. Laure, Y. Demay, and J. 
Polym. Process. 19, 118 (2004). 
R. Valette, R. Laure, Y. Demay, and J. 
Polym. Process. 18, 171 (2003). 

C. D. Han and R. Shetty, Pol. Eng. Sci. 18, 180 (1978). 
R. Valette, P. Laure, Y. Demay, and J. F. Agassant, Jour- 
nal of Non-Newtonian Fluid Mechanics 121, 41 (2004). 
K. L. Helton and P. Yager, Lab. Chip 7, 1581 (2007). 
A. Utada, E. Lorenceau, D. Link, P. Kaplan, H. A. Stone, 
and D. Weitz, Science 308, 537 (2005). 

P. Guillot, A. Colin, A. S. Utada, and A. Ajdari, Phys. 
Rev. Lett. 99, 104502 (2007). 

D. J. Korteweg, Arch. Neerl. Sci. Ex. Nat. II (1901). 

R. B. Bird, C. F. Curtiss, R. C. Armstrong, and O. Has- 

sager, Dynamics of polymeric liquids, 2nd ed., Vol. 1. 

Fluid Mechanics (Wiley, New York, 1987). 

P. Petitjeans, C. R. Acad. Sci. 322, 673 (1996). 

C. E. Hickox, Phys. Fluids 14, 251 (1971). 

M. d'Olce, J. Martin, N. Rakotomalala, D. Salin, and 

L. Talon, J. Fluid Mech. 618, 305 (2009). 

K. Chen and D. D. Joseph, J. Non-Newton. Fluid Mech. 

42, 189 (1992). 

Y. Amarouchene, D. Bonn, J. Meunier, and H. Kellay, 
Phys. Rev. Lett. 86, 3558 (Apr 2001). 



